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Abstract 

We present simulation results for lattice QCD with light pions. For 
the quark fields we apply chirally symmetric lattice Dirac operators, in 
particular the overlap hypercube operator, along with the standard over- 
lap operator for comparison. This allows us to simulate at very low pion 
masses. The results are related to Random Matrix Theory and to Chiral 
Perturbation Theory in order to extract information about the pion decay 
constant, the scalar condensate and the topological susceptibility. 

1 Chiral Perturbation Theory 

When a continuous, global symmetry breaks spontaneously, we obtain 
a continuous set of degenerate vacuum states. Expanding around one 
selected vacuum, one distinguishes between excitations to higher energy 
(which are identified with massive particles) and fluctuations, which pre- 
serve the ground state energy. The subgroups of the energy conserving 
symmetry group can either transfer the selected vacuum to a different 
vacuum state, or leave it simply invariant. The number of generators re- 
lating different vacuum states corresponds — according to the Goldstone 
Theorem — to the number of massless Nambu- Goldstone bosons (NGB) 
that emerge. At low energy, the NGB can be described by an effective 
theory as fields in the coset space of the spontaneous symmetry break- 
ing (SSB). Such effective theories still apply if we add a small explicit 
symmetry breaking; we then deal with light quasi-NGB, which dominate 
the low energy physics. The effective Lagrangian C c g contains terms of 
the quasi-NGB fields, which obey the original symmetry, as well as the 
(explicit) symmetry breaking terms. All these terms are hierarchically or- 
dered according to some low energy counting rules for the momenta and 
the quasi-NGB masses. 

This concept is very general, but it was introduced in the framework of 
chiral symmetry breaking in QCD. At zero quark masses QCD is assumed 
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to exhibit a chiral SSB of the form 

SU{N f ) L ® SU(N f ) R -» SU(N f ) L+R , (1) 

where Nf is the number of flavours involved. In this case, the coset space 
is again SU(Nf), and the corresponding low energy effective theory is 
known as Chiral Perturbation Theory (xPT) PP- 

A small quark mass supplements a slight explicit symmetry breaking, 
and the quasi-NGB are then identified with the light mesons, i.e. the pions 
for Nf = 2 — and for Nf — 3 also the kaons and eta particles. 

In view of our lattice study, we have to put the system into a finite 
volume; we choose its shape as V = L 3 x T (T > L). Then xPT was for- 
mulated in two regimes, with different counting rules for the terms in £ e ff- 
The usual case is characterised by Lm, 3> 1, where is the pion mass, 
i.e. the lightest mass involved, which corresponds to the inverse correla- 
tion length. This is the p-regime, where finite size effects are suppressed, 
and one expands in the meson momenta and masses (p-expansion) 

The opposite situation, Lm, < 1, is denoted as the e-regime. In that 
setting, an expansion in the meson momenta is not straightforward, due to 
the important role of the zero modes. Fortunately, the functional integral 
over these modes can be performed by means of collective variables P|. 
There is a large gap to the higher modes, which can then be expanded 
again, along with the meson masses (e-expansion) 011]. 

In both regimes, the leading order of the effective Lagrangian (in Eu- 
clidean space) reads 

£efi = ^Tr[9 M (7 t 9 M [/]-isTr[X([/ + (7 t )] + ... , 

U £ SU(Nf) , M = dia S (m u ,m d , (m s )) . (2) 

The coefficients to these terms are the Low Energy Constants (LEC), and 
we recognise F and E as the leading LEC. As it stands, E and F occur 
in the chiral limit M = 0; at realistic light quark masses F turns into 
the pion decay constant F w . Experimentally its value was measured as 
F w ~ 93 MeV, which is somewhat above the chiral value of F ~ 86 MeV 
[7j. E is not directly accessible in experiments, but its value is assumed 
to be >(250 MeV) 3 . 

The LEC are of physical importance, but they enter the xPT as free 
parameters. For a theoretical prediction of their values one has to return 
to the fundamental theory, which is QCD in this case. There is a notorious 
lack of analytic tools for QCD at low energy, hence the evaluation of the 
LEC is a challenge for lattice simulations. 

The LEC in nature correspond to their values at V = oo, and the 
p-regime is close to this situation. However, it is interesting that these 
infinite volume values of the LEC can also be determined in the e-regime, 
in spite of the strong finite size effects. Actually one makes use exactly 
of these finite size effects to extract the physical LEC. Generally, we need 
a long Compton wave length for the pions, l/m n , and in view of lattice 
simulations in the p-regime we have to use an even much larger box length 
L. In this respect, it looks very attractive to work in the e-regime instead, 
where we can get away with a small volume. 

However, such simulations face conceptual problems: first, light pions 
can only be realised if the lattice fermion formulation keeps track of the 
chiral symmetry. In addition, the e-regime has the peculiarity that the 
topology plays an important role 0: X?T predictions for expectation 
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values often refer to distinct topological sectors, so it would be a drastic 
loss of information to sum them up. 

2 Lattice QCD with Chiral Fermions 

These conceptual problems could be overcome only in the recent years. 
The solution is the use of a lattice Dirac operator D which obeys the 
Ginsparg- Wilson relation Its simplest form reads (in lattice units) 

D75 + 75,D= -D 75 D , n>l, (3) 

which means that D _1 anti-commutes with 75, up to a local term that 
vanishes in the continuum limit. Even at finite lattice spacing, this local 
term (275 /n in eq. (J3J) does not shift the poles in the propagator D _1 . 
This fermion formulation has a lattice modified, but exact chiral symmetry 
[H] and exact zero modes with a definite chirality. That property provides 
a definition of the topological charge by means of the fermionic index 
v = n+ — n_ (where n± is the number of zero modes with positive/negative 
chirality) 0. 

The simplest solution to this relation is obtained by inserting the 
Wilson-Dirac operator Dw into the so-called overlap formula |10|. 

D%} = A ov + n , A ov = fiAo ( AlAo J , A = D - n ■ (4) 

H. Neuberger suggested this solution with Do = Dw, and we denote 
the resulting D y as the Neuberger operator (at mass zero). We used 
it in all our applications presented below at fi = 1.6. It is motivated, 
however, to study also the generalisation with different kernels Do, in 
particular when they already represent an approximate Ginsparg- Wilson 
operator We suggested to use a kernel with couplings in a unit 

hypercube on the lattice (hypercube fermion, HF) , which are constructed 
with Grassmannian block variable renormalisation group transformations 
|12|. Its gauging also involves "fat links" |13| (the same is also true for 
the alternative HF kernel of Ref. |14p. Compared to Dw, multiplications 
with this kernel require more numerical work (about a factor of 15 in 
QCD), but part of it is gained back immediately since the evaluation of 
the overlap operator — which has to be approximated by polynomials in 
practical implementations — has a faster convergence. There remains an 
overhead of about a factor 3, but there are further gains of the resulting 
overlap HF in terms of locality, rotation symmetry and scaling. These 
virtues have all been tested and confirmed for free fermions |11) and for 
the 2-flavour Schwinger model |15|. In QCD we worked out such HF 
kernels at (3 = 6/ go = 6 1131 and recently also at f3 = 5.85 |16| . which 
corresponds to lattice spacings of a ~ 0.093 fm resp. a ~ 0.123 fm (in 
quenched simulations). Again an improved locality — see Fig. — and 
rotation symmetry could be confirmed, while a systematic scaling test 
is still outstanding. The axial anomaly is correctly reproduced in the 
continuum limit of any topological sector for the Neuberger fermion |18|. 
as well as the overlap HF |19|. 

A number of subtle numerical tools for simulations with overlap fer- 
mions have been elaborated |20) . but their simulation is still computation- 
ally expensive. For the time being, only quenched QCD simulations are 
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Figure 1: The function f(r) measures the maximal correlation between sink and 
source if they are separated by a distance r. The decay of its expectation value 
must be (at least) exponential for locality to hold. Our plots show that this is the 
case for both, the Neuberger fermion \17J] and the overlap HF, at f3 = 6 (on the 
left) and at (3 = 5.85 (on the right) on L 4 lattices. We also see that the decay 
is clearly faster for the overlap HF, i. e. its locality is significantly improved. 



possible; the inclusion of dynamical quarks 21 might be a challenge for 
the next generation of supercomputers. 

In view of the evaluation of Low Energy Constants, the quenched 
approximation causes logarithmic finite size effects |22|. Nevertheless it is 
important to explore the potential of this method, in particular in view 
of the evaluation of LEC in the e-regime. 



3 The Pion Mass 

For Wilson's traditional lattice fermion formulation, it was a insurmount- 
able problem to reach light pion masses; due to the additive mass renor- 
malisation and other conceptual problems, the pion masses always re- 
mained above about 600 MeV (at least quenched). In order to verify how 
close we can get to realistic pion masses with the overlap HF, we first 
evaluate the corresponding pion mass in the p-regime, as a function of 
the bare quark mass m q (which we assume to be the same for all flavours 
involved). The latter is added to the chiral operator D&t of eq. (|1J as 




The evaluation considers the exponential decay of the pseudoscalar cor- 
relation function, which is deformed to a cosh behaviour by the periodic 
boundary conditions. The method can still be improved by subtracting 
the scalar correlations function, i.e. we study the decays of both, 

C P (t) = ^(P t (x,t)P(0)) , and 
C PS (t) = ^<P f (x,t)P(0) - S\2,t)S(0)) , (6) 

X 

where P and S are the pseudoscalar and the scalar density. In Cps some 
contamination by topological finite size effects is eliminated. The result is 
shown in Fig. (on the left) , which illustrates that the data are in agree- 
ment with the expected behaviour mjj. cx m q . At m q = 0.01 (in lattice 



units) we arrive at a value of ~ 230 MeV, i.e. significantly closer to the 
physical value than the simulations with Wilson fermions. At this point, 
we are already close to the transition to the e-regime, m^L m 1.7. The 
corresponding pion masses for the Neuberger fermion as measured by our 
collaboration with more data can be found in Ref. |23|. A measurement 
in exactly the same way (using the same 29 configurations) based on Cp 
is shown for comparison in Fig. [5] on the right. 
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Figure 2: Left: the pion mass vs. the bare quark mass (both in lattice units) for 
the overlap HF, evaluated from the temporal decay of the pseudoscalar correlator, 
with or without subtraction of the scalar correlator. The lowest pion mass in 
this plot (at m q = 0.01) corresponds to about 230 MeV. This is a preliminary 
result, based on 29 propagators. On the right we show the comparison to the 
Neuberger fermion, both evaluated from Cp in exactly the same way. 



4 The Topological Susceptibility 

In the quenched approximation, the topological charge (which is identified 
with the fermion index v) does not depend on the quark mass. Its statis- 
tical distribution is expected to be Gaussian. Fig. |^] (on the left) shows 
part of our index history for the overlap HF and the Neuberger fermion. 
They deviate a little, (|^ v-hf — z/n|) ~ 0.89(4). The plot on the right 
shows, as an example, our index histogram measured with the Neuberger 
operator at /i = 1.6, which is reasonably consistent with a Gaussian. 

From the point of view of the naive, constituent quark model, the 77' 
meson is amazingly heavy (m v i ~ 958 MeV). An explanation for this 
property was given based on topological windings of the Yang-Mills gauge 
field. In particular, the Witten-Veneziano formula relates m n i to the 
topological susceptibility Xt in the leading order of a 1/N C resp. Nf/N c 
expansion (N c is the number of colours), 

"V = "j^f* 4 ' X* = 77 > v '■ topological charge. (7) 

In this formula \t i s understood as a property of pure gauge theory 
(hence it is sensible to define it with the quenched fermion index) , whereas 
the other terms refer to full QCD. A recent large-scale study with Neu- 
berger fermions |24| arrived in the continuum extrapolation at \t = 
(191 ± 5 MeV) 4 , or — in dimensionless units — Xt r o = 0.059(3) (where 
rrj is the Sommer scale, which translates lattice quantities into physical 
units). This is consistent with a heavy r\ meson. 
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Figure 3: On the left we show part of the index histories for the overlap HF 
and the Neuberger fermion, for the same configurations on a 12 3 x 24 lattice at 
(3 = 5.85. The plot on the right shows the distribution of topological charges on 
a 16 3 x 32 lattice at f3 = 6, measured with the index of the Neuberger operator 
at fj, = 1.6 on 506 configurations. We see a decent agreement with a Gaussian 
distribution. Its width represents the topological susceptibility, see Fig. ^[ 



In Fig.0]we show our preliminary result for Xt r o m a physical volume of 
(1.48 fm) 3 x 2.96 fm at lattice spacings a = 0.123 fm and a = 0.093 fm. We 
include results for the Neuberger operator, as well as the overlap HF, and 
we also show the continuum extrapolation of Ref. |24) for comparison. Our 
present data suggest a trend to a somewhat larger value of Xt r 0i though 
the overlap HF is closer to the result of the literature. More statistics will 
be required to clarify if we arrive at a result consistent with Ref. [24 | . 
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Figure 4: The topological susceptibility Xt r o °^ different lattice spacings in the 
same volume V = (1.48 fm) 3 x 2.96 fm. On the left we show our result with the 
standard statistical evaluation. Alternatively we cut the measured distributions 
at \v\ = 1.5, 2. 5,..., 5. 5 and matched them to a Gaussian each time, so that 
in each case a larger statistics is involved. The results obtained in that way is 
shown on the right. So far our value for Xt is a little higher than the one of Ref. 
\24\ , which is, however, based on a larger statistics. 



6 



5 Determination of E 



Chiral Random Matrix Theory (RMT) introduces a pseudo Dirac operator 
/ iW \ 

Drmt = I w | n I , W : complex n+ x n_ random matrix. (8) 



v iW t 

In this way, QCD is simplified to a Gaussian distribution of fermion matrix 
elements. The hope is to capture nevertheless some properties of QCD 
in the e-regime. In fact, the Leutwyler-Smilga spectral sum rules [Sj were 
successfully reproduced in this way |25|. 

There are also explicit conjectures for the probability distribution p 
of the low lying eigenvalues i\ of the continuum Dirac operator. One 
introduces the dimensionless variable z = AEV and considers the spectral 
density 



V P ( :Hz) . (9) 



v— — oo 



Formally this term is taken at V = oo (so one deals with a continuous 
spectrum), although the prediction refers to the e-regime. Note, however, 
that the quark mass is set to zero. This shows that the conjecture involves 
a number of assumptions, and a test against lattice data is motivated. The 
last step in eq. © is a decomposition into the contributions of the topo- 
logical sectors, which only depend on \u\. In a fixed sector we can further 
consider the individual densities of the leading non-zero eigenvalues, 

n=l,2,3,... 

Chiral RMT provides explicit predictions for the leading densities pn (z) 
pa), which are supposed to hold up to some energy (depending on the 
volume) . These predictions have been compared to lattice data for overlap 
fermions in Refs. |27) . which did agree for not too small volumes, up to 
some value of 2. 1 The only free parameter involved in these fits is the 
scalar condensate E. The successful fits provide therefore also a value for 
E, which is found in the range where it is generally expected. At least 
(2) agreed well with the RMT prediction at \v\ = 0, 1, 2 for L > 1.1 fm. 

In Fig. we show a new result, for the mean value of the leading 
non-zero eigenvalue, (Ai), 2 in the topological sectors \v\ = 0, 1, . . . , 5 on 
different lattices with different overlap operators, but again in a fixed phys- 
ical volume. We then fit the RMT predictions for the optimal value of E to 
these data, which yields values in the range E = (268 MeV) 3 . . . (290 MeV) 3 . 

At this point we mention that also staggered fermions have a remnant 
chiral symmetry on the lattice, though not with the full number of gener- 
ators (in contrast to Ginsparg- Wilson fermions). Thus they are also free 
of additive mass renormalisation. The standard formulation turned out 
to be "topology blind", but a sensitivity to \v\ — similar as the one shown 
in Fig. — can be obtained by a suppression of the mixing between its 
pseudo-flavours |28) . 



^The work by Giusti et al. is the most extensive study in this context. 

2 The spectrum of our Ginsparg- Wilson operators is actually located on a circle in the 
complex plane with radius and centre fi. For a comparison with the RMT predictions in the 
continuum, we map this circle with a Mobius transform onto the imaginary axis. 



7 



0.05 
0.045 

0.04 
0.035 
a 0.03 
v 0.025 

0.02 
0.015 

0.01 
0.005 

1 2 3 4 5 

I V | 

Figure 5: We show our data for the leading non-zero Dirac eigenvalue (Ai) 
(resp. (zi) — T,V(Xi)) in the topological sectors \v\ = 0, 1,...,5, in a volume 
(1.48 fm) 3 x 2.96 fm. For the Neuberger operator we measured at lattice spacings 
a = 0.093 fm and a — 0.123 fm, and in the latter case we also used the overlap 
HF operator. The optimal fit leads to a value for £ in each case, which is 
indicated in the plot. 

6 Mesonic Correlation Functions 

To obtain direct predictions for quenched simulations results with chiral 
lattice fermions, also quenched has been worked out for mesonic 

correlation functions [29]. The vector correlator vanishes to all orders, 
and the scalar and pseudoscalar correlators involve additional LEC in the 
leading order, which are due to quenching. This is not the case for the 
axial vector correlator, where a parabolic (rather than a cosh) behaviour 
was predicted in the e-regime. In fact, this shape could be observed 30 
if the spatial box length exceeds about L>1.1 fm, i.e. the same bound 
that we encountered in Section 5. This also allowed for a determination 
of F-X. On the other hand, the sensitivity of the curve to E is too weak 
to evaluate it in this way. We further noticed that the Monte Carlo 
histories in the sector v — are plagued by strong spikes, which are related 
to the significant density of very small Dirac eigenvalues |30| (at higher 
charges, this density is suppressed). Hence a brute force measurement in 
the topologically neutral sector would require a tremendous statistics (see 
also Ref. |31|). A strategy to avoid this problem was proposed and applied 
in Ref. |32j . Work is also in progress to establish a lattice gauge action, 
which preserves the topological charge over long periods in the Monte 
Carlo history 1331. This would help us to measure expectation values in 
a specific sector — which is desirable in the e-regime. With the standard 
plaquette gauge action (that we are using so far) it is tedious to collect 
statistics in a specific topology. 

Here we turn our attention to another procedure, which only considers 
the zero-mode contributions to the correlators. Following Ref. |34| we focus 
on the pseudoscalar correlator, where some re-definitions allow us to study 
only F and one quenching specific LEC called a in the leading order of 
quenched x?T. We distinguish the connected and the disconnected zero- 
mode contributions, 

C H (x) = {v){x)vi{x)v\ (0)«j(0)) , V H (x) = (vUx)vi(x)vt(p)vi(p)) , ( n ) 
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Table 1: The statistics of configurations at \i/\ = 1 and 2, for different lat- 
tices and overlap operators, which was included in our study of the zero-mode 
contributions to the pseudoscalar correlator. The result is shown in Fig. 



where i,j are summed over all the zero-modes, D^vt — 0. It is difficult 
to fit these quantities directly, but it is easier to fit the data to the leading 
order in the expansion around the minimum in the time coordinate t [34] 
(after summation over the spatial lattice sites x). This minimum is at 
T/2, hence we expand in s — t — T/2, 



L*ds {1 



T 1 



(12) 



where we nee lect C((s/T) 3 ). The combined fit of our data in the sectors 
\v = 1, 2 to C and T> — at different values of s — leads to the results 
for F and a shown in Fig. IS1 fll)]. see also Ref. |35|. (In these fits we also 
made use of our measured results for {v 2 ).) Since eqs. (1 1 2> only hold up 
to 0((s/T) 3 ) we should extrapolate down to small s. We only obtained a 
stable extrapolation, with acceptable errors, for the overlap HF, but not 
for the Neuberger operator on two different lattice spacings — although 
the statistics is similar, see Tabled 

The result of these extrapolations agrees within the (sizable) errors 
with Ref. |34| . Note that these are values for the bare parameters. In 
particular the renormalisation of F involves the factor Za, which amounts 
to about 1.45 in the chiral limit for our parameters of the Neuberger 
operator at /3 = 5.85 |23| . 
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Figure 6: Results for F and a (LEC in the leading order of quenched \PT), 
based on our data for the zero-mode contribution to the pseudo-scalar correlation 
function in the topological sectors \v\ = 1 and 2. In both cases, we obtain a stable 
behaviour and a decent extrapolation down to a small fitting range s only for the 
overlap HF. The statistics for these plots is given in Table QJ 
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7 Conclusions 



We discussed the potential of relating simulation results for lattice QCD 
with chiral fermions to the predictions by x?T. For the first time we 
applied the overlap HF in this context and compared its results also with 
Neuberger's standard overlap operator. We gave preliminary results for 
the pion mass in the p-regime. In the e-regime we are interested in the 
determination of the LEC that appear in the effective Lagrangian. 

This project is still on-going, as part of the efforts of the xLF col- 
laboration to simulate QCD close to the chiral limit. For the approach 
discussed here we saw that we can go down close to 200 MeV with the 
pion mass. The evaluation of the LEC in the leading order seems feasible, 
though a larger statistics is required for this purpose. So far we can report 
that data obtained in the e-regime — i.e. in small volumes — can indeed 
be matched with the analytical predictions by chiral RMT and quenched 
XPT. We expect the same property also for dynamical chiral quarks, once 
our machines are powerful enough to simulate them. 
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